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Abstract 

We show how to locate all the transmission eigenvalues for a one 
dimensional constant index of refraction on an interval. 



1 Introduction 

The scattering operator for the time harmonic Helmholtz equation with com- 
pactly supported index of refraction n(x) maps the asymptotics of incident 
waves (Herglotz wave functions) v to the asymptotics of (outgoing) scattered 
waves u + . Both u + and v° are defined in M. n and satisfy 

(A + k 2 n 2 ) u + = (1 - n 2 )k 2 v 
(A + k 2 ) v = 

in all of R n . 

If the scattering operator has a zero eigenvalue, we say that k 2 is a trans- 
mission eigenvalue [B], [I]. It follows from Rellich's theorem and unique 
continuation that the scattered wave u + is identically zero outside any do- 
main D that contains the support of n, and therefore that u = u + and v = v 
are a pair of nontrivial functions satisfying 

(A + k 2 n 2 )u = (1 -n 2 )k 2 v in D (1.1) 

(A + k 2 )v = in D (1.2) 



u 



= p 

dD du 



dD 
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Whenever such a nontrivial pair exists, we say that k 2 is an interior 
transmission eigenvalue for the pair (D,n)Jj 

Remark 1. If we scale v by introducing w = k 2 v , then the k 2 on the right 
hand side of equation U.l\) disappears, and we see that the interior transmis- 
sion eigenvalue problem is a true generalized eigenvalue problem which can 
be written as 

Au — (1 — n 2 )w = —k 2 n 2 u 
Aw = —k 2 w 

with the same boundary conditions as in U.3\) JMj. 

The analogous interior problem for obstacle scattering is the Dirichlet 
problem, 

(A + k 2 ) v = in D 
v =0 

dD 

and the Dirichlet eigenvalues are fairly well understood. In particular, the 
Dirichlet eigenvalues of a one dimensional interval of length L are exactly 

2 2 

^7?-. The one dimensional interior transmission eigenvalue problem, even 
for a constant index of refraction, is not as simple. In one dimension, 
with the index of refraction n(x) equal to the positive constant a, equations 
(TO) fTO|) fTQj) become 

u" + k 2 a 2 u = k 2 (l-a 2 )v (1.4) 
v " + k 2 v = (1.5) 
ii(±j) = ii'(±f) = (1.6) 

so that k 2 is an interior transmission eigenvalue if and only there exist a 
nontrivial pair of eigenfunctions (u,v) satisfying f ll.4pfll.5pfll.6p . 

We will show in proposition H]that k 2 is an interior transmission eigenvalue 
if and only if A; is a root of the equation 

(a - I) 2 sin 2 ((a + _ {a + if sin ^ a - 1)|) = (1.7) 



*See section 8 of [3], as well as [8], for more details about the connection between 
transmission eigenvalues and interior transmission eigenvalues 
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and the two theorems below describe the roots of (11.71) . which we will call the 
interior transmission wavenumbers (positive square roots of interior transmis- 
sion eigenvalues). 

Theorem 2. The interior transmission wavenumbers satisfy 

|T .... 2, 3cr + l 

lm(k < -log, 

L \cr — 1] 

In theorem [3] below, the notation [^J denotes the floor, the greatest in- 
teger less than or equal to — . The condition 1-^-1 = 1 1 means that there 

° ^ m LmJ L m J 

is no integer between and — , and Z denotes the integers. 

° mm" ° 



Theorem 3. Let m = and let j e Z 

i. If j- ^7* and Z, i/iere are exactly two simple interior transmis- 

sion wavenumbers in the strip 

^ <Re(ib)<^ + 1 ^ 



(a+l)L w (cr+l)L 

If L^-J = L^"J ^ e ?/ are complex conjugates with nonzero imagi- 
nary parts 

2. If — G Z, t/iere a quadruple interior transmission wavenumber at 
k — — and no others in the strip 

m 1 



2(j 



1 )i/ 



1)1/ 



cr 



The transmission wavenumbers for two cr's, with L = 2, are plotted below. 
The vertical blue lines indicate the sets Re(fc) = c^^[\ L ■ The roots that fall 
on these lines are quadruple. 
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As these cr's are rational, the transmission eigenvalues are periodic. Re- 
mark |5] below explains how we did the computation. 
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2 Lemmas and Proofs 

Because we are dealing with constant coefficients, the transmission wavenum- 
bers k are the zeros of the analytic function on the left hand side of equation 
(11.71) . This analytic function has two factors, and our proofs simplify if we 
treat each factor separately. 

Proposition 4. k is a transmission wavenumber if, and only if, k satisfies 
(a - 1) sin((a + 1)^) = [a + 1) sm((a - 1)^) (2.1) 

or 

(a - 1) sin((a + 1)|) = -(a + 1) sin((a - 1)^) (2.2) 

The wavenumber k satisfies Ii2.1\) if there exists an eigenfunction pair (u, v) 
of odd functions of x, and k satisfies 112. ty) if there exists an eigenfunction 
pair (u, v) of even functions. 

Proof. If a pair u, v satisfies (II .4111 .6p . then so do their odd and even parts, 
and the calculation is a little simpler for each of these parts. The odd parts 
satisfy (ll.4til.5l) with (II. 6p replaced by 

M (0) = v (0) = u(f)=« / (f) = 
and the even parts satisfy (ll.4til.5l) with (II. 6p replaced by 

u '(0) = 1/(0) = «(f ) = u'(f ) = 
In the odd case, the boundary conditions at zero imply that 

v = sin(kx) 

and 

u = — sin(fcx) + Asm(kax) (2.3) 
for some constant A. The boundary conditions at x = ~ imply that 



sin(crA;— ) cos(£;— ) = a sin(&— ) cos(aA;— ) (2.4) 
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Noting that the left and right hand sides are the beats representation of the 
combination of the frequencies (cr+1)^ and (a— l))|f suggests that we rewrite 
flUD as 

(a - 1) sin((a + = (a + 1) sm((a - 1)-*) 

which is (12. ip . A similar calculation shows the existence of even transmission 
eigenfunction pairs if and only if k satisfies (12. 2p . □ 

If we introduce the new variables 

m 

and 

z = (a 

then equations (12.1 ft and (12. 2\\ become equations for a z whcih depend on the 
single real parameter m. They become 

p° m {z) := sin(m7T2;) — msin(7T2;) = (2-6) 
Pm( z ) := sin(m7rz) +msin(7T2;) = (2.7) 

where the superscripts refer to even and odd, respectively. Because a > 0, 
\m\ > 1. Because equations (12. 6ft and (1 2. 7ft remain the same if we change the 
sign of m, we may also assume that m > 1, and will do so throughout)]]. 

17TZ 

Remark 5. If m — 2 is rational and we set w = e <? , then ([OP becomes a 
polynomial equation in w. 

w 2 P-l = l ( w p+i - yjP-i) (2.8) 
q 

We can use a polynomial root finding algorithm to find the solutions and then 
take logarithms and rescale to find the interior transmission wavenumbers k. 
This is how we produced figure 1 . 

With the notation introduced in (12.51) . theorem [2] becomes 

' Alternatively, we could start with the assumption that a > 1, and then use the 
relationship between the transmission wavenumbers of a and those of kj(^) = akj(a) 
to show that our formulas and estimates apply to a < 1 as well. 



a+1 



a 



l)k 



L 

2^ 



(2.5) 
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Theorem 6. If z satisfies ( TO]) or (27fy , th 



en 



, , , log(2m + 1) , 

Im(z <-2j 2.9 

7r(m — 1) 

Proof. \i z = x + iy satisfies (12. 6p 

e m " l?/l - e" m " l2/l < 2| sinmvr^l = 2m| shi7rz| < me l|sl + me" l|!/l 
so that 

e mn\y\ < e -m7r|y| + me ^|y| + n^-Tlwl 

Dividing both sides by e 71 "' 2 ^ gives 

e (m-i)ir\y\ < e -(nH-i>|»l + m + me -2^\y\ < 2 m + 1 

Now taking the logarithm of both sides yields ( 12.9ft □ 

Locating the real parts of the roots of p m (z) = Pm( z )Pm( z ) wm take more 
work. We will show (roughly) that there are always two roots in certain of 
strips of the complex plane. We define 



I± = < x 

m 

Sj = 



j ^J + l 

— < x < 

m m 

— < Re(z) < 

m m 

A restatement of theorem [3l using the notation of equation (I2.5p . is the- 
orem [7| below. 

Theorem 7. Let m > 1. 

1. If neither 3- nor ^±1 are integers, then p m has two simple roots in the 
strip S± . 

m 

(a) If the interval I± contains an integer k, both roots are real, and 

m 

one is bigger and the other smaller than k. 

(b ) If the interval Ij_ does not contain an integer, the roots are com- 

m 

plex conjugates (with nonzero imaginary parts). 

2. If 3- is an integer, — is a quadruple root of p m , and p m has no other 
roots in S j or Sj-i . 



Theorem [7] is an immediate consequence of the analogous theorems for p° m 
and p e m . These theorems are a little more complicated to state, but easier 
to prove. 

Theorem 8. Let m > 1. 

1. Suppose that neither i- nor ^±1 are inteqers. 

rr mm J 

(a) If the interval I± contains an integer k and j + |_— \is even, then 
p° m has a real root in the interval (—,k) and no other roots in S± 
and a root in (k, and no other roots in S±. If the interval 
Ij_ contains an integer k and j + \_— J is odd, then p° m has a real 
root in the interval (k, ^—) and no other roots in S j and pt, a 
root in ( — ,k) and no other roots in S j . 

(b ) If the interval I± does not contain an integer and j + |_— J is odd, 
p° m has 2 simple complex conjugate roots in S± and p e m has none. 

m 

If j + L^J is even, p e m has 2 simple complex conjugate roots in 
S j_and p° m has none. 

2. If — is an integer, neither p° m nor p e m have any roots in S± or Sj^i . If 

m mm 

j + [^\ is even, — is a triple root of p° m and a simple root of p e m . If 
j + L^J is odd, — is a triple root of p e m and a simple root of p° m . 

We will prove theorem|H](we will only prove the results for p£j by studying 
a one parameter family of equations with m fixed and parameter (5. 

Pp(z) = sm(rmiz) — /3sm(iiz) (2-10) 

We will use the fact that that the roots of the entire function pp(z) de- 
pend continuously on the parameter We will also use the fact that, as 
j3 — > oo, the roots of p@ must approach those of sin(7T2;). This last fact 
follows from the first by multiplying f)2.10p by a = 4, and the fact that the 
roots of a sin(m7T2;) — sin(7rz) depend continuously on a as a approaches zero. 



*The continuity of the roots of analytic functions follows from the fact that roots are 
isolated, together with the argument principle. 
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Before giving the rigorous details, we give a rough outline of the proof of 
Theorem [HI To simplify this outline we ignore the case where — is an integer. 

Theorem [8] makes two assertions 

1. It tells us the number of roots p° m in each strip S±. This number is the 

m 

same for pp for all j3 > 0. 

2. It tells us the number of real roots in each strip when [3 = m. The 
number of real roots of pp changes as (3 increases from zero to m, but 
is the same for all f3 > m . 

We establish the number of roots in each strip by observing that: 

1. at /3 — all roots are real and on the boundary of the strips S±, and 

m 

the derivatives are nonzero and pointing in the direction of the correct 
strip (lemma [9]), so that, for small positive /3, the correct number of 
roots have entered each strip (corollary [TO]) . 

2. Lemma [LTJ tells us that no roots can cross the boundary of Sj_, so the 

m 

correct number of roots remain in each strip for all (3 (corollary [T2l). 
We calculate the number of roots which are real at f3 = m by observing that: 

1. As (3 approaches infinity, a single root must approach each integer (the 
roots of sin %z) , and the other roots must leave every compact subset 
of C. 

2. A calculation(lemma [T3|) shows that for f3 > m, all roots are simple. 

3. This means that any real roots of pp in S± at (3 = m must remain real 

m 

for (3 > m and remain in S±, so they remain in the interval I±. But, 

m m 

as j3 approaches infinity, the roots of p@ must either approach the roots 
of sin7rz or leave every compact subset of C. The roots in the S±'s 

m 

that don't contain an integer must leave every compact set, so they 
cannot remain in the interval I±, so they cannot have been real when 

m 

(3 = m. 

4. The S±'s that contain an integer contain only one root for all j3, so the 
root remains simple and therefore real at (3 = m. That root approaches 
the integer as (3 approaches infinity. 
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The detailed proof follows: 

Lemma 9. The roots ofpo = sm(m7iz) are —. Each of these roots is simple, 
and thus continues to a unique root Zj(/3) of pp for small (3. Moreover, 



dZj_ 

d/3 



= (_iy+L^J| sin (^[)| (2.11) 
,3=0 m 



Proof. At f3 = 0, the roots of smmnz are located at — and are all simple. 
The implicit function theorem implies that, for small /?, they remain simple 
and depend analytically on /3, so we can define Zj(f3) to be the continuation 
of that root. Differentiating ( 12.10)) yields 

(mir cos(mnZj) + (3n cos^nzj)) = sm(irzj) 

so that 

■IV" sin & = f— lV+L^Jl sinfiZL)! 



d8 



p—0 cos it j v / m v ' 1 * m < 



□ 



Knowledge of the derivatives tells us that the roots are in the correct 
strips for small f3. 

Corollary 10. For /3 > and sufficiently small, all roots of pp are real and: 

1. Suppose that neither j- nor 2±1 are inteqers. 

rr mm 3 

(a) If the interval Ij_ contains an integer k, then pp has a real root in 
the interval {^,k) and no other roots in Sj_ if j + [^\is even, pp 
has a root in (k, ^-) and no other roots in S± if j + [^\is odd. 

(b ) If the interval I± does not contain an integer and j + |_— J is odd, 

m m 

Pp has no roots in S± . If j + |_— J is even, pp has two roots in Sj_ . 

m m m 

2. If 3- is an integer, j^is a simple root of pp. If j + is even, there is 
a simple root in each of I± and h-\ . If j + [— J is odd, there are no 

mm m 

roots in I ±_ or I ,-i . 
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Proof. Because sinz = sinz, the roots of ( 12. 6 p must occur in conjugate pairs, 
and simple roots must remain real for small j3. Lemma M tells us that, Zj 
starts at — and moves to the right if — is not an integer and 7 + 1 — 1 is 
even, or to the left if — is not an integer and 7+1 — 1 is odd. If I j does not 
contain an integer, then + I 2±i I j s oc [d whenever 7 + 1^-1 is even, and 
even whenever j + |_— J is odd, so that both Zj and Zj + \ move into I± when 
j + |_— J is even, and away from I± when j + )_— J is odd, which establishes 
itemllbl 

If I± contains an integer k, then Zj moves into Ij_ and Zj+i moves away 

m m 

from Ij_ when j + [— J is even, and Zj \ moves into I± and Zj moves away 
when 7+1 — 1 is odd. The root which moves into / j remains in either ( — ,k) 
or (k, for small (3, which establishes item [Taj 



Finally, if — is an integer, both terms in (I2.10p vanish, so that ^ remains 
a root for all f3, and remains simple at least for j3 small. If — is an integer, 
t± + lirlj = j + I ± \ _ 2 and *±! + I i±l| = j + \ ±\ + l^so moves 

m L m -I J LmJ m LmJ •' LmJ ' J 1 

into f j-i and Zj+i moves into J j+i if j + [— J is even, and both move away 

from those intervals if j + |_— J is odd. This establishes item [2] and finishes 
the proof. □ 

At a root of pp(z), the real parts of the sinm7r,2 must equal the real part 
of /3 sin nz. The lemma below identifies the sets on which the real parts of 
these terms vanish. 

Lemma 11. Re(sin(z)) = iff sin(Re(z)) = 
Proof. 

sin(x + iy) = sin x cosh y + i cos x sinh y 

and cosh y never vanishes. □ 

A consequence of lemma [TT] is that pp has no roots on Ke(z) = — or on 
Ke(z) = k, when j or k are integers, unless — equals an integer. This is 
enough to show that the roots which entered each strip for small /3, must 
remain there. 
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Corollary 12. For all < f3 < oo 



1. The conclusions in items [7al andUM of corollary [Wj remain true. 

2. If — is an integer and j + |_^J is odd, — is a simple root of pp and pp 
has no other roots in Sj_or Sj+i . If — is an integer and j + |_— J is 
even, pp has three roots in Sj-i U S±, one of which is —. 

mm m 

Proof. The hypothesis of item [TJ that — and ^— are not integers, combine 
with lemma [TTJ to guarantee that pp can have no roots on the boundary of 
Sj_. As the roots of the entire analytic function pp depend continuously on 

m 

j3, this implies that no root can enter or leave S±. Thus, if S± contained no 

m m 

roots for small (3, it contains no roots for all (3. If S± contained two roots 
for some f3, it contains two roots for all f3. This establishes item [lb] for all /3. 
To establish item [Taj note that lemma [TTJ forces the simple root to remain in 
k < Re(z) < or — < Re(z) < k. As the simple real root remains simple, 
it must also remain real. 



The two statements in item [2] follow similarly from the fact that no root 
can cross — or — , as these cannot be integers if — is an integer because 

mm' ° m ° 

m > 1. □ 

Lemma 13. For (3 > m, all roots of pp are simple. For f3 = m, all roots are 
simple except z = — in the case that j + — is an even integer. In this case,— 
is a triple root of pp. 

Proof. All roots of pp that aren't real are simple, because there are no more 
than three in any strip, and they must occur in conjugate pairs. Suppose we 
have a double real root, then 

sinm7T2; = f3smnz (2-12) 

and 

mcosmixz = f3cosnz (2-13) 
Squaring both f)2.12p and f)2.13p . and adding the results gives 
sin 2 rrniz + cos 2 muz + (m 2 — 1) cos muz = f3 2 

or 

cos 2 m^=£^ (2.14) 
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which is impossible for (5 > m because the right hand side is bigger than one. 
If (5 = m, then ( I2.14p and (I2.13P are possible, but only if 



cos7r^ = cosmirz = ±1 

which implies that z and mz are either both even or both odd integers. This 
is equivalent to the statement that z = — and j + [^\ is an even integer. 
One can check directly that the root is indeed triple in this case. □ 

Proof of Theorem^ We only prove the theorem for ]f ml as the proof for yf m 
is analogous. Item [la] follows immediately from item [1] of corollary [121 Item 
[1] of corollary [12] also implies all of item [Tbl except the fact that the two 
roots are not real. To see this, note that, as j3 — > oo, the roots of pp must ei- 
ther approach the roots of sin irz, which are integers, or leave every compact 
subset of the complex plane. Since J x contains no integers, the two roots in 

m 

Si must leave every compact subset of C. As Ke(z) must remain between 

m 

i- and ^— , the imaginary parts of these two roots must grow unbounded 
as (3 — > oo. If they were real at /3 — m, lemma [13] would imply that they 
were simple and remained simple and real for all (3, contradicting the fact 
that their imaginary parts grow unbounded as /3 — > oo. Hence they must not 
have been real at = m. 

Item [2] follows for a similar reason. We can check directly that ^ is either 
a simple or a triple root, depending on the parity of j + |_— J- Item [2] of 
corollary [12] tells us that these are all of the roots in the closure of Sj-i US±, 

m m 

so there are no other roots of p^in these strips. □ 

3 Some Limiting Cases and Conclusions 

We have shown that there are two transmission wavenumbers in each closed 
strip (Jj^ L < Re(fc) < t^jy^ of width j^jjj;, and two real transmission 
wavenumbers in each of those strips which contains an integer. We show 
below that, in the limiting cases of a weak scatterer (a — > 1), there are 
no longer any real transmission wavenumbers, and in the limiting case of a 
strong scatterer (a — > oo), all the transmission wavenumbers become real. 

The Born, or linear approximation, is the limit as a — > 1. We return to 
equation f ll.7p . set nz = kL, and let a — > 1 to obtain 
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ttz = i sin nz (3-1) 

Introducing the parameter j3 in front of the ttz on the right hand side, 
and repeating the steps we used to prove theorem [3, gives: 

Theorem 14. For every integer j except j = —1 and j = 0, equations A3. 1\) 
have exactly two non-real complex conjugate roots in the strip j < Re(j') < 
and they lie in the interior of the strip. The only root in —1 < Re(z) < 1 
is the triple root at z = 0. 

The imaginary parts of these roots grow logarithmically. The precise 
asymptotics, which can be verified by substituting (13.21) into (13. ip . are 

4 = 0' + §)±* log(2|i + ||) + O (3.2) 
so that kj = ^ satisfy 

kf~i ((j + i)±<log(2|i + |D) 

where, for the odd integers j, Zj solve equation (13. ip with the plus sign, and, 
for even integers j, the Zj solve (13. ip with the minus sign. 

The limit as a — > oo is a little more work to calculate. This is the limit 
as m — > 1, so that equation (12.71) becomes 

sin nz = — sin ttz 

or 

sin7T2; = (3.3) 

with roots at the integers. Thus, as a — > oo, the even transmission wavenum- 
bers 

3 aL 

have the same asymptotics as the square roots of the eigenvalues of the 
Dirichlet problem 

u" + k 2 a 2 u = 
u(±§) = (3.4) 

13 



with odd (u(x) = —u(—x)) eigenfunctions. 



The solutions of (12. 6p . which describe the transmission wavenumbers with 
odd eigenfunctions, have a different asymptotic behavior, that will match the 
asymptotics of a different self-adjoint boundary value problem. To see the 
limiting equation, we subtract sin nz from both sides of (12. 6|) 

sin(m7rz) — sin nz = m sin nz — sin nz 
and divide by m — 1 

sin(m7T2;) — sinnz 



m — 1 

which becomes, as m — > 1 

nzcosnz = sinnz 

or 

nz = tannz 
The transmission wavenumbers satisfy 



m sin ttz 



k(a - 1)^ = tanfc(<r - 1)^ (3.5) 

and thus have the same asymptotics, as a — > oo, as the square roots of the 
eigenvalues of the self-adjoint boundary value problem. 

v" + k 2 a 2 v = (3.6) 
±f t/(±f ) = V (±f ) 

with odd eigenfunctiong§. 

Because these limiting wavenumbers are the square roots of the eigen- 
values of a self-adjoint problem, we see that, as a — > oo, all transmission 
wavenumbers become real. 

In our discussion of the one dimensional transmission eigenvalue prob- 
lem, we have observed certain relationships between transmission eigenvalues, 

§The wavenumbers of (13. 6[) with even eigenfunctions satisfy equation ()3.5|) with tangent 
replaced by cotangent 
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Dirichlet eigenvalues, and the eigenvalues of another self-adjoint boundary 
value problem. We expect (optimistically) that some of these relationships 
will persist in higher dimensions. We know that, in the Born approximation, 
there are no real transmission eigenvalues [8], that a generic radially sym- 
metric index of refraction has complex transmission eigenvalues [7], and that 
there are infinitely many transmission eigenvalues [S] [Tj . All of these results 
are consistent with our one dimensional conclusions. 



A consequence of theorem [3] is that the counting function for the number 
of interior transmission wavenumbers in a strip of width W is 

7T 

which has the same asymptotics as the counting function for Dirichlet wavenum- 
bers for an interval of the same length with index of refraction n — a + 1. 

The counting function for real interior transmission wavenumbers in an 
interval of width W is 

7T 

which matches the asymptotics of the counting function for Dirichlet wavenum- 
bers for an interval of the same length with index of refraction n = a — 1@ 
It seems reasonable to conjecture that the relationship between the counting 
functions for the Dirichlet wavenumbers and the transmission wavenumbers 
remains true for non-constant a, and in higher dimensions as well. 
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